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УДК 539.1 М-16 

ѴАПП!Т1Н'': РЕЗОНАНС в мамторазбавленных твердых телах при низ¬ 
ких ТЕМПЕРАТУРАХ: Препринт ИТЭФ 87-10/ 

Ф.С.Джепаров, Е.К.Хеннѳр - М., ШІИИатоминформ, 1987 - 23 с. 

Рассмотрен парамагнитный резонанс в неупорядоченной систе¬ 
ме дшіольно-взаимодѳйствугащих спинов При произвольных темпера¬ 
турах зѳемановской подсистемы и при условии, что где 

$ - обратная дипольная температура, а Ё-энѳргия взаимодейст¬ 
вия спинов на среднем расстоянии. Сформулированы кинетические 
уравнения, вычислены теплоемкости подсистем, а также в рамках 
модели Андерсона получено явное выражение для функции формы ли¬ 
нии. Показано, что при условии ^ Е 0 ^>1, где Е 0 -энѳргия взаимо¬ 
действия ближайших спинов, термодинамика спиновой системы, в от¬ 
личие от высокотемпературного случая, в основном, определяется 
взаимодействиями на средних расстояниях. Для расчета формы ли¬ 
нии предложен метод корреляционного разложения, а для расчета 
теплоемкостей использовано разложение по концентрации. 

У.З.ВяНерагоѵ, Е.К.Неппег 

МАОНЕТІС КЕ50НАІІСБ ІН МАОЯЕТІСАЬЬУ ПУШТЕВ ЗОІДЦЗ АТ ЬО'Я ТНЕИРЕ- 
КАТОНЕЗ. 

ТЬе рагатайпеііе геаопапев іп іНе йізогПвгвй вріп вуеіетв 
ѵчііЬ аіроі-йіроі іпівгасііопв ів сопвіавгвд. ТЬе іНеогу 1в ехіеп- 
в^а іпіо іЬе гееіоп, ѵНеге ^ Е4У (^?іа іЬе іпѵвгвв ііроі іЬвтрв- 
гаіиге, Е ів іЬе іпіегасііоп епегву аі іЬе теап аіаіапсе). ТЬв 
кіпвііс ециаііопа агв Гогтиіаіеа, ала іНе Ьеаі сарасіііеа Гог гваі 
вузіета апа іііе гевопапсв Гогш Гипеііоп Гог іНе Апавгвоп тоаві агв 
саісиіаіеа. Іі ів эНочт іЬаі іЬв ѳріп №е гто йупаті с в аі ^3 Е 0 
(Е іа іЬе іпіегасііоп впвг^У Гог іНе пвагеаі пеііНЪога) іа ѵегу 
аіГГегепі Ггош іНе Ні^Ь іЬетрегаіиге опе апа іа гийеа Ъу іЬе іп- 
іегаеѣіопа аі аѵагаве аізіапсев. ТЬв соггеіаііоп еірапаіоп ів 
ааѵапсеа іо саісиіаіе іНе Гогш Гипоііоп апа іЬв Ьваі сарасііівв 
аге саісиіаіеа иаіп§ іЬе сопсепігаііоп ехрапэіоі.. 

Рис. - , список лит. - 21 найм. 
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.),1Ѳ87 



I. ВВЕДЕНИЕ 


Изучение величин, характеризующих магнитный резонанс в сис¬ 
теме дипольно взаимодействующих спинов, является одной из самых 
фундаментальных задач радиоспектроскопии. С необходимостью рас¬ 
чёта или, хотя бы, оценки формы линии, термодинамических величин 
и скорости магнитной релаксации сталкиваются практически все тео¬ 
рии, направленные на изучение свойств многочастичных систем. В 
случае магниторазбавленных веществ для оценки ширины и формы ли¬ 
нии использовались, в основном,два подхода: метод моментов и ста¬ 
тистический метод Андерсона [і,2] . Первый был подвергнут серьёз¬ 
ной критике [3,4] , а второй, после ряда обобщений [б-?) , при¬ 
годен для описания поглощения энергии высокочастотного магнитно¬ 
го поля слабонеравновесной спиновой системой в высокотемпературном 
приближении (ВШ) , когда Т 00 ° С 60 * -центральная час¬ 

тота магнитного резонанса, Сі) ц -частота локального поля, обус¬ 
ловленного диполь-дипольными взаимодействиями). Термодинамические 
свойства таких систем и магнитная релаксация в ВШ изучены весь¬ 
ма детально Г 1,2] . 

Для последнего времени, однако, характерен интерес к изуче¬ 
нию спиновых систем в области низких температур, вплоть до 
магнитного упорядочения [в-І0| . При этом используется динамичес¬ 
кое охлаждение резервуара спин-спиновых взаимодействий,причем 
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в соответствии с теорией Провоторова спиновая система находится в 
состоянии двухтемпературяого квазиравновесия. Указанная теория, 
однако, высокотемпературна по отношению как к зеемановской под¬ 
системе, так и к подсистеме взаимодействий, и построение низко¬ 
температурной теории насыщения упирается в соответствующие тео¬ 
рии формы линии, магнитной релаксации и термодинамических свойств. 

Распространение статистической теории формы линии на низко¬ 
температурный случай было проведено для некоторой частной ситуа¬ 
ции в работе [іі} в работе [і2] изучалось влияние спин-поляриза- 
ционных эффектов на фазовую релаксацию. Подсистема взаимодейст¬ 
вий в этих работах подразумевалась высокотемпературной. Шаги в 
сторону построения низкотемпературной теории насыщения были сде¬ 
ланы в работах [іЗ-Іб] , однако ни в одной из них не проводился 
анализ формы линии и термодинамических величин ни для регулярных, 
ни для разбавленных систем. Изложенный же в Цв] анализ нелинейных 
эффектов в низкотемпературной термодинамике налагает столь силь¬ 
ное ^граничение на температуру, что для неупорядоченных систем 
он практически бесполезен. 

В данной работе развит метод корреляционного разложения, кото¬ 
рый позволил построить статистическую теорию формы линии в магни¬ 
торазбавленной системе при произвольном охлаждении зеемановской 
подсистемы (т.ѳ. высокой спиновой поляризации) и на той стадии 
охлаждения дипольной подсистемы, для которой /ЗЕ« I, где 1 - 
энергия взаимодействия спинов на среднем расстоянии, а уЗ -обрат¬ 
ная дипольная температура. С помощью концентрационного разложе¬ 
ния в том же приближении вычислены теплоёмкости указанных подсис¬ 
тем, входящие в уравнения теории насыщения. 
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В конкретных расчетах мы ограничивается случаем, когда спин 
5* =1/2, образец имеет эллипсоидальную форму, а распределение 
спинов по узлам решетки нѳкоррелировано. 

2. СИГНАЛ ПОГЛОЩЕНИЯ ПРИ ДВУХТЕШЕРАТУРНОМ 

КВАЗИРАВНОВЕСИИ. УРАВНЕНИЯ ТЕОРИИ НАСЫЩЕНИЯ 
Гамильтониан рассматриваемой системы имеет вид 


#= 76 л + 'К. я ~со 0 $ 2 , Сі) 

где -секулярная часть гамильтониана диполь-дипольных взаи¬ 

модействий, й> 4 = С (Ц^-амплитуда вращающегося вокруг оси % 
поля частоты со ). В условиях квазиравновесия матрица плотности 

3- (~Н 'Х л )/Вр*хр(-р Я '-ряз*). (2) 


Общее выражение для сигнала поглощения в этих условиях на рас¬ 
стоянии А от ларморовской частоты в приближении линейного 
ответа представимо в форме И 


(д)> — - е < [&+ (4), •$.]> 


( 3 ) 


ще временная 


определяется оператором 


Для дальнейшего существенен переход от коммутаторного сред¬ 
него С «)=< [3+ &.*.]> к антикоммутаторному /І(і)= <ІЗ+(^),3.] + У , 
Представим (2) в виде $-=:*хр( г р9С)/&рехр(-р9ё) .где 
96 • 

Перейдем в 3+ Сі) к гамильтониану : 

5+ (і) = &хр(~; ё^ео с -б) $Е-6)5+&хрСг;Э*-і) ■ 



После этого шеем: 


<4 (*» = - с&)* С *~ (4) 

где С(6) формально соответствует равновесию с температурой /в 
и гамильтонианом Рі. . Используя известную связь спектральных 
представлений равновесных корреляционных функций С СО и 
А (*) [16] и переходя после этого от РС к исходному гамиль- 

тониану^ получим 

<$,(*)>= , 

*Й>= І&жі) _5* а<«~рО»*) ■ 

Эти соотношения обобщают флуктуащшнно-диссипациояную теорему 
на двухтемпературяый случай. Переход от (3) к (5) важен тем, 
что функция д (л) неотрицательна при всех Л . Это доказы¬ 
вается записью д(д) в собственном представлении гамильтониана 
. Таким образом, сигнал поглощения при произвольных 
температурах подсистем обращается в ноль в единственной точке 
4 = ё^.со 0 Св ВТП это свойство хорошо известно). Приближённая 
процедура вычисления < 3^ (А)> не должна нарушать этого 
важнейшего свойства. Из него, в частности, следует вожсжность 
измерения при любом соотношении р и путём реги¬ 

страции нуля сигнала поглощения. 

Выражение (5 ) мажет стать основой для получения уравне¬ 
ний теории насыщения магнитного резонанса в твердых телах, обоб¬ 
щающих уравнения Провоторова, справедливые в ВШ, на случай 
произвольных значений р и р х . Используя точное соотношение 
[8]. вычисляя <^ а > с помощью С2) 
и считая, что временная эволюция системы описывается величинами 
/2 Х (4) и р (-6), получим с помощью (5) 




где С л * 77 ? < \*^ ^-пардиаяъные теплоёмкости подсистем 

ал * ъ р* 

( ^ • 

В качестве второго уравнения, как и в ВТП /8/, можно при¬ 
менить закон сохранения энергии спиновой системы во вращавшейся 
системе координат. При этом гамильтониан 2€ вск =А , а 

матрица плотности^выраженная через зешановскую температуру в 
лабораторной системе, по-прежнему имеет вид (2) (дипольная тем¬ 
пература в обеих системах одинакова). Условие ^ ~>- О 


( < -‘Ы + Ѣ 0 С 2сі)~^ + С^х 4 ' и)с С гг) ^ - О • ( 7 ) 

Из (6) г. (7) находим 

^ Ц.ЫВіА- Я-V*$.(&), 

<и~ х*/ СялС&-С**& л *- 


& = т - (Ѳ) 

х+/ С хх С лл - с гл > 

Мы ограничились здесь сравнительно элементарным выводом 
этих уравнений. Отметим однако, что при более последовательном 
анализе на основе проекционной техники Накаджима-Цванцига с 
использованием проектора Кавасаки-Гантона \Т7] получается тот 
же результат. 

Уравнения (8) в ВТП переходят в уравнения Провоторова. 

Они отличаются от уравнений, сформулированных в [Іб} отсутст¬ 
вием неоправданного, на наш взгляд, переопределения подсистем. 



В отличие от ВТП и $(&) являются функциями и /3, 

которые в общем случае найти в явном виде невозможно. Ниже с 
помощью корреляционного разложения найдено д (&) , а с помощью 
концентрационного разложения - . 

3. КОРРЕЛЯЦИОННОЕ РАЗЛОЖЕНИЕ СИГНАЛА ПОГЛОЩЕНИЯ 


Важным допущением статистической теории формы линии являет¬ 
ся возможность ограничиться лишь анизотропной частью гамильтониа¬ 
на 3^ (модель Андерсона): 

^ " *^ п ***Рі* > 4к =/ Ъ*Фк , е$в’ 3елХ %д,АиГ°- (9) 

Здесь п- -число заполнения узла ^ : П- = і, если узел } 

занят спином, и П- = 0, если не занят. 

Объектом ближайшего рассмотрения является корреляционная 
функция /460= <&[.$+ (*), 5_] + . Здесь означает 

конфигурационное усреднение по всем возможным реализациям распре¬ 
деления спинов по узлам решетки. Используя операторное тождество 


-л* д * «Ч*- 5 ** 




в справедливости которого можно убедиться, например, сравнивая 
матричные элементы в 3^ - представлении спиновой пары, при¬ 
ведём /)&) К виду /](*)= 11<п к <П&ср (і/і- к 5 *^). 

Непосредственное вычисление А (^ ) не представляется воз¬ 
можным. В данной работе развит метод, намеченный ранее в [з] и 
состоящий в последовательном учете многочастичных корреляций, 
аналогичный в известной мере групповым разложениям теории газов. 

В рассматриваемой задаче можно выделить два типа корреляций: 
динамические и статистические. Специфика модели Андерсона - воз¬ 
можность точного учета динамических корреляций. 
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Перепишем А (*) переходя от среднего по матрице плотности 
(2) к средним по "зеемановской" матрице плотности 
8с - °хр(гА*. ХгУЗр ѳхр <**.), 

М±>^<п к <;*ф[1 г п.(іА }к $* -Д*-$-у_п п А ;т )]> в / (іі) 

/< ахр[-А 21 п. (А; к -+і;І_*п т Л] т )]^ X . 

В (II) явно выделены те члены взаимодействия , которые 

обуславливают парные корреляции спина " К ". Слагаемые же, 
входящие в 21 , обуславливают корреляции более высокого поряд¬ 
ка, и поэтому разложим А ( ^ ) по числам заполнения . 
Разложение функции чисел заполнения по этим числам 

имеет вид 

Ф({п})= Ф м +Лп к (4 с %)- Ф") + 

+ іТ^п^Ф ( %,<*у4> ( %<уФ Со) У — * (І2 ) 

- Ф 1 % Л Л п к п к --п к Л С-*) Ф 

/П=І П, ‘ <,—<т ' п р-о 

Здесь Ф <0) соответствует всем незанятым узлам, - одному за¬ 
нятому и т.д. Эти формулы можно получить, например, используя 
тождество <?(•>)+ п^СО-ФС 0 )) и метод математической 

индукции. В сумме по <,•••>>, следует исключить любые совпадения 
значений переменных суммирования и учесть, что все суммирования 
ограничены одним общим объемом. Перед переходом к бесконечному 
объему выражение под знаком указанной суммы надо симмѳтризовать. 

После конфигурационного усреднения разложение по числам 
заполнения приводит к разложению искомых макровеличин по степе¬ 
ням концентрации занятых узлом <^~ і*'- число спинов). 



и при достаточно ограничиться первыми членами разложения. 

Применительно к (II) разложение по п т является корреляционным 
разложением и в первом порядке даёт 
й(*)• А А %, А 21 <п к А*%)> с , 

АцТ/. п Г<^&А&&~ М<МФо -А М МЪ (ІЗ) 

* * Рс ’ 


/ 4 * <*ф[Т.»;(іАх2*4 А^)}>сА м й(-А^ г !/А;х)> € , * (і4) 

На каждом шаге корреляционного разложения термодинамические сред¬ 
ние берутся по сравнительно небольшому числу спинов и могут быть 
в принципе найдены точно, хотя объем вычислений быстро растет с 
ростом порядка разложения. Отдельную проблему, которую также 
удается решить, составляет конфигурационное усреднение. 

Вычисления квантовостатистических средних при 3* =1/2 мож¬ 
но гшолнить на основе хорошо известных свойств 6~ - матриц. 

В частности, имеет место тождество 

г«Г> . ш 

Преобразуем с его помощью числитель в (14): 

х”( і) - <«*>Гр!;<Ч< *?*-м>***а )]>.= . 

= <[р<е^(;»■/!;, ф -р> л і » 

=5' л+р (ч^сч-аЗШ >•= (І6) 


+ 0*-р, і-+-і) ■ 


9 

Здесь зеемановское" среднее по состоянию спина ^ , 

р_ - спиновая поляризация при бесконечной 

дипольной температуре. Обозначая знаменатель в (14) через %'с , 
имеем X™ (+=°) . Таким образом, А { %)='Е.<п к Х С е У^ < с>^ 

Используя далее легко проверяемое замещением Пу на 1 или 
О тождества 

П-ЛЬЛ , о4>$іг>$=. і+Пу(оЬ^-і) , (17) 

ввделим в X? и общий множитель П[4+/1;(сб(!(&-4)2, 
сократим его и получим 

Г<л* > е , 

X^Й> = Ц? 1р1*гГ>;(4- СЛ; К (і-рЩ к )-1+іЛ; к (р~Щк))р 

+ Гр- -р, с— 4 -; , * ». х . {І8) 

Здесь ал. к = ап = 4$ . 

Используя тождество 

П^і-Пу (А;**;)]- (71(4+ ПуАу) (р(4+ , (19) 

ввделим в А (ві (і) множитель 

&К (*)= О &&&) = П^у (<Ы> М ~І+ 'РЩ*)], ( 20 ) 

обеспечивающий точный учет динамических корреляций, подвергнем 
оставшееся выражение дополнительному разложению по числам 
заполнения и получим основной и первый поправочный члены: 

У'йу= /> пю ( ф >* /)"%)- Г<*. - 

К 



( 21 ) 



ІО 


Выражение (21) допускает конфигурационное усреднение. При некорре¬ 
лированном распределении спинов по узлам <п- 

и т.д. После усреднения имеем 

А 10,г> Ь)-= 4(<^Ьк~і^Р 4*0*)] , 

х у+е (22) 

А (0 ^С+) = -ір(і-р х )4 Х Т- П ' 

Кт? /?Кт 

и.- *>х~ __ . С23) 

і~р Х ^кт 

Выражения (22), (23) в определенном смысле являются окончатель¬ 
ными. В них вычислены как термодинамические, так и конфигурацион¬ 
ные средние. Далее, для получения точных количественных результа¬ 
тов, следует конкретизировать форму образца, тип кристаллической 
решетки, ориентацию внешнего поля относительно образца и произво¬ 
дить численные расчеты по формулам (22), (23). Этот путь, однако, 
малопригоден как в практическом, так и в теоретическом плане. 

Оценки показывают, что нельзя рассчитывать на непосредственное 
вычисление полученных решеточных сумм для сколько-нибудь "большо¬ 
го" образца за приемлемое время. В теоретическом отношении из 
сравнения (22) и (23) нельзя непосредственно установить истин¬ 
ный параметр, по которому произведено дополнительное разложение; 
маловероятно, что таким параметром является . Как показывают 
многие прецеденты, ответ на эти качественные вопросы и выполне¬ 
ние асимптотических оценок выражений типа (22), (23) можно полу¬ 
чить в приближении сплошной среды (ПСС), к которому мы перейдем 
после того, как сменим порядок (^) по ^. 

Вычисление термодинамических средних, входящих в А ^60, 
производится с помощью тождества 

&], ( 24 ) 



II 


получаемого двухкратным использованием формулы С15). Преобразова¬ 
ния, аналогичные произведенным выше при вычислении (^), 
позволяют привести ( ^ ) к виду 






(25) 


где 


+ (26) 
-о--Л ■ 

Здесь обозначения сл,- х , , ^;к имеют тот же 

смысл, что и в (18), а = />4^4)». Выделяя в (25) общий 
множитель (4) и дополнительно разлагая оставшееся выраже¬ 
ние по числам заполнения, получаем, что А^ ~ 

Таким образом, (23) полностью содержит члены порядка 
в разложении (13). 

4. ПРИБШИШИЕ СПЛОШНОЙ СРЕДЫ ( ПСС ). 

Из формул (22) и (23) видно, что ПСС для функции формы 
линии справедливо при выполнении условий /«I и Е 0 -6 » "І , 
где . а -наименьшее расстояние между 

спинами в решетке. Первое из этих условий позволяет пренебречь 
запретом на совпадения координат разных спинов, а благодаря 
второму условию в той области, где функции сильно меняются при 
переходе из узла в узел (и поэтому, вообще говоря, замена сумм 
на интегралы незаконна) они одновременно сильно осциллируют и 
вклад этих осцилляций зануляется. Ограничимся рассмотрением 
эллипсоидальных образцов в высокосишѳтричных (например,кубических) 
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решетках. Для них фигурирующие ниже решеточные суммы вида 
У практически не зависят от положения узла х 

Переход к ПСС в (22) осуществляется следующим образом: 

П (4+/•(;*)= еоср X- 

где С = -объемная концентрация спинов ( П - число узлов 

в единице объема). Центр интегрирования соответствует узлу "к " 
Полный переход в (22), (23) к ПСС дает 

А (о,о \±)^Ме с0±С *\ (<гр х )сМй х (*)* С<Хі ^І (28) 

где 




і- р *им@ 


. />т* 


Вычисление интегралов вынесено в Приложение I. Используя (П4), 
(П5), имеем ,е, 

/)« 0 , С+)=Уе. *' Ъ)%і 

(в длинновременном пределе). Здесь 2> — Лж х с^ х ^/СЗ'Гз)у 
$=[^(/Г-4,&2)-7г4]рсг І '4 ; ('-/<§<4 В зависимости 

от ориентации эллипсоида относительно внешнего поля и соотноше¬ 
ния длин осей). Видно, что параметром, по которому ведется кор¬ 
реляционное разложение во временном представлении^является ве¬ 
личина % -~/&Е , т.ѳ. произведение обратной 
температуры подсистемы взаимодействий на энергию взаимодействия 


спинов на среднем расстоянии. 



13 


Как показано в Приложении I, в Фурьѳ-гірѳобразовавиѳ величи¬ 
ны д(0Д) (-6 ) вносят большой вклад относительно малые времена. 
Поэтому для получения сигнала поглощения уместно выполнить 


Фурьѳ-прѳобразованиѳ непосредственно от А' 


,С0,І) 


( 4 ), определен¬ 


ной формулами (28) - (30). Используя (П8), в делом для функции 
фор^ получаем 

(ЗІ) 

Отметим, что /XI < І/2.Простейшее условие применимости (31) 
можно сформулировать как 


/«/«*,/^ 4"*^. (зіа) 

Однако, формула (31) остается физически осмысленной при всех 
/Ѳ а и А . Поэтому возможно, что первое из неравенств (Зіа) 
слишком сильное и в действительности оно излишне. Для более 
определенного вывода необходим анализ следующего члена корреляци¬ 
онного разложения. 


5. КОНЦЕНТРАЦИОННОЕ РАЗЛОЖЕНИЕ ТЕРМОДИНАМИЧЕСКИХ 
ВЕЛИЧИН 

Энергии и теплоемкости можно выразить через функцию 

Я&А)=- <А 5р е~ Л = **> л > с , 

где /^= Р(о,р л ) . Очевидно, что 

<ж+ >= у с *«- • 
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Разлитая р по числам заполнения из и производя усред¬ 
нение, в главном порядке по имеем 

Р ш Е* ~ ^ ^к* = Р л +г с , ( 32 ) 

Здесь 96^ СЗ, к) есть гамильтониан дипольного взаимодействия 
двух спинов, расположенных в у' и к . Запишем 96^ О,К) 


* •$. 4 ). (ад 

Стандартному дипояь-дипольному гамильтониану отвечает Д = 3 -, 
а модели Андерсона 9 =0 . Учитывав, что » X 3 = &■> 

применяя лето проверяемую формулу 

*<■%?«> р(І)+ $•?(- ?> Х1&&У ГС- Ю], -2- і . 

и соотношение (24), можно получить, что 


<е 




Вообще говоря, расчет теплоемкостей (кроме С хх - (і-^) ) 
весьма громоздок и необходим численный анализ. Отметим лишь, “ 
что и С гх не зависят, а С^ х зависит от формы образца. 
Существенные упрощения происходят в высокотемпературном прибли¬ 
жении, когда р Е 0 <& I ив противоположном, "низкотемпера¬ 
турном" случае, когда »-І. Напомним, что при этом 

по-прежнему ~^Е е <ъг 1 . 

Результаты ВТП легко получаются разложением полученных вы¬ 
ражений по ^ .В "низкотемпературном" случае как Г 
так и энергии 'Э Р/Э/З л не допускают перехода к пределу 
сплошной среды, т.ѳ. замены сумм на интегралы по всему прост¬ 
ранству, вследствие сингулярности А- к на малых расстояниях. 
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Однако, для теплоемкостей такой переход возможен и приводит к 




> с Мг = Я V* Г(Р) - 


(35) 


При этом для V . зависящей только от Р, получается представ¬ 
ление однократным интегралом (подробности см. в Приложении 2). 

Эти результаты означают, что в температурной области, где 
Е теплоемкости,в основном, определяются взаимодей¬ 

ствиями на средних (а не малых, как в ВТП) расстояниях, и, соответ¬ 
ственно, термодинамика оказывается совершенно иной, чем в ВТП. 

6 . ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ 

Рассмотрим выражение для формы линии (31). Как уже отмеча¬ 
лось, и, следовательно, д. (Л)~>0 , что необходимо 

согласно п.2. В случае высокотемпературной подсистемы взаимодей¬ 
ствий, когда для любых представляющих интерес частот 
т.ѳ. при , (31) і 


$(А)~ 


(36) 


уг + 

В этом случае можно сопоставить наши результаты с результатами 
работ [іІ,І2]. Как и там, мы приходам к лоренциану, ширина кото¬ 
рого Л) не зависит от поляризации. Однако, даже для сферического 
образца формула (36) предсказывает существование сдвига частоты 
магнитного резонанса, пропорционального концентрации и поляри¬ 
зации спинов. Этот эффект, который на первый взгляд противоречит 


классической макроскопической электродинамике, целиком опреде¬ 
ляется флуктуациями в случайном распределении спинов; обсуждение 
этого вопроса проведено в работе [ів] . 

В случае, когда система полностью размагничена, а дипольная 
подсистема охлаждена. 



іе 

< *“ ^ ( 36 ) 

т.е. сигнал поглощения почти антисимметричен относительно лармо¬ 
ровской частоты (с точностью до сдвига) и обращается в ней в ноль. 
Подобное поведение сигнала поглощения хорошо известно в случае 
высокотемпературной подсистемы взаимодействий [8] . 

Полученные формулы допускают (как и в ВШ) использование 
отрицательных дипольных температур.,Это свойство верно в любом 
конечном порядке концентрационного разложения и обусловлено 
ограниченностью спектра 76^ для любого конечного числа спинов. 
Однако, в отличие от ВШ, вообще говоря, ъ/&Уф- Г/ в). 

То обстоятельство, что параметром малости, соответствующим 
корреляционному разложению, оказались величины х. и % , делает 
выражение (31) применимым к изучению широкого круга явлений,свя¬ 
занных с динамическим охлаждением подсистемы взаимодействий. 

В частности, оказывается доступной температурная область, в кото¬ 
рой начинается переход системы дшольно-взаимодействующих 
спинов в фазу дипольного стекла [іэ] . 

В следующем порядке корреляционного разложения весьма гро¬ 
моздкие расчеты, связанные, в частности, с вычислением двухцѳнт- 
ровкх интегралов, показывает, что дяиняовремешая асимптотика 
(-Лз/ёІ+Сѣ*)) и поэтому в низкочастотной 
асимптотике ширина линии Л и сдвиг 5" получают приращения, 
пропорциональные ^.Простая картина высокочастотной асимптотики 
нам неизвестна. 

Интересно отметить, что корреляционное разложение в любом 
порядке свободно от объемных расходимостей, появление которых 
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столь усложняет работу с выеокотемдературиши разложениями 
теории магнетизма [8] . Это можно показать в общем виде исполь¬ 
зуя (і2). 

Как отмечалось вше, модель Андерсона позволила произвести 
точный учет динамических корреляций для функции формы линии. 
Аналогичная теория в МП дает вполне удовлетворительные резуль¬ 
таты. Однако выход за пределы ВТП в связи с охлаждением резер¬ 
вуара взаимодействий ставит вопрос о роли изотропной части га¬ 
мильтониана, включающей флип-флопы, при низких температурах. 
Развитый в данной работе метод не распространяется непосредствен¬ 
но на этот случай. Формально это связано с невозможностью факто¬ 
ризовать ,<ц. (ё) . Простое концентрационное раалечениѳ /}(ё) , 
возможное для произвольного взаимодействия, напрямую непригодно, 
т.к. в каждом его порядке возникает незатухающая функция времени. 
Разработан, однако, другой подход к теории функции формы линии 
в разбавленных системах [2СІ] , базирующийся на том же разложении 
по числам заполнения (12) с привлечением метода функций памяти. 
Этот подход позволил решить задачу о форме линии в ВТП с учетом 
изотропной части взаимодействия. Распространение этого метода на 
случай низких температур позволяет сделать предварительное заклю¬ 
чение, что при ± отличия от результатов, получен¬ 

ных в модели Андерсона, невелики. Этот вопрос впоследствии будет 
обсужден подробно. 

Найденные в данной работе выражения для формы линии и тепло¬ 
емкостей достаточны для решения уравнений (б) , описывающих 
кинетику насыщения спин-системы в твердых телах, в указанном 
вше температурном диапазоне. Однако, не менее важные аспекты 
теории насыщения связаны с учетом сяин-решѳточной релаксации. 
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В отличие от ВТП простой феноменологический ее учет был бы неубе¬ 
дителен, т.к. "времена релаксации" вне ВТП зависят от температур 
подсистем. Этот вопрос требует самостоятельного изучения. 

Другой очень важный вопрос состоит в том, что в широком 
диапазоне амплитуд переменного поля неупорядоченные систе¬ 

мы ведут себя как неоднородноуширешгае. Это выявлено в работе 
/21/ и должно приниматься во внимание при сопоставлении изложен¬ 
ной выше теории с экспериментом. 


Приложение I 

Рассмотрим вначале 0^ (4) . Представим его в виде 

, где 


*(■&) » ({- Л* , Т(*)=&\^г> 


А(?)± _ 


Имеем: /Къ) = ЪГ^Сі~2сг^)/СІ ъ 3 ). Обозначим 

сделаем замену переменных Х=І/г, 3 , У= с»9 и получим 

"■ ? . т 

Здесь Х^ - , А ЛЛЛ = . Для тех интегралов, кото¬ 
рые рѳгуляризованы на малых расстояниях, -» «=» , 

для бесконечного сферического образца Х^ и . После интег¬ 

рирования по частям и использования асимптотических представлений 
интегрального синуса и косинуса получим, что для сферического 
образца при ІП^ Х //=%* і 

ѵі/= , 7"= (ПЗ) 

ЗГЗ О 


Етодящий в (ПЗ) интеграл равен но О [ІІ^, а ~~ (/X- <5и ). 
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Таким образом, в этом приближении й±(Ь) = (з,80І4І~ 0, 5~0 2 р-6) . 

В случае конечного образца эллипсоидальной формы сумма, порожда¬ 
ющая интеграл И/ , 

О 

составлена из быстро сходящихся решеточных сумм 21 
и от местонахождения узла К практически не зависит. 

Решеточная суша, порождающая интеграл Т, 

содержит вклад от первого слагаемого в правой части, зависящий, 
вообще говоря, от формы образца и местонахождения узла к ; 
для эллипсоида, однако, 2ГА Х от местонахождения "< " не за¬ 
висит. В этом случае представим 

г ип фГ = 21(«, , 

и первая 21 перестает зависеть от формы образца и может быть 
вычислена по сфере с центром в " к ". Учтем, что для кубических 
решеток ^8_3 -Лтгігк, где -і<: 5^2. (в зависимости от 

соотношения осей эллипсоида и ориентации внешнего поля; для сфе¬ 
ры =0). Таким образом, для эллипсоида 

с %&)= ЯМ-&-6, СП4) 

где Ъ = ^сг-%%=[^(Г?-&^)-7Гъ]рсг г * . 

Оценим . Ограничимся вначале простой ситуацией высо¬ 

котемпературной подсистемы взаимодействий. Тогда 

4Й)* (П5) 
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Отказ от условия высокотѳмпературности не изменит длинновремен- 
ной асимптотики 2 Л (6) , определяемой (П5). Для оценки доста¬ 

точно ограничиться аппроксимацией 4И а. = а при /а/< I и 
41п а = -А. при /а/> X. Переходя к переменным х,у и разби¬ 
вая область интегрирования в плоскости Ху на части,соответствуйте 
"высокотемпературным" и "низкотемпературным" спинам,получим, что 
вклад в от "низкотемпературных" спинов пропорционален 
и > т * ѳ * П Р И йольших временах, -б-*»/*, , 

мал по сравнению с (Д5). 

это не означает, что вклад от "низкотемпературных" 
спинов мал в фурьѳ-прѳобразовании величины , которое 

интересует нас в связи 


это Фурье-прео Сразованиѳ непосредственно: 

А^(л)~ А С6 ‘^(+) е. --;р €. 


;**-*%&) 


ЪС+), (Пб) 


где О х Оі) определяется формулой (30), а О^С-Ь) -формулой (П4). 
ИМѲШ 


А‘ А %)=рО-р*)еУ№ 


^-(з/)(ъ)-Д~Ъ) я -+2> Х 

Переходя к координатам ■2, Э, <А и делая указанные при 

выводе (П2) замена переменных, получим, выполнив интегрирование 

по У, 


і-рЧА м&р 

^_7 




^ ■ тг 7&Г 


№В) 
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где /4= -^-(а+%), -Р- % ЛЭ. Предст а в им подынтегральную 

функцию в виде ^&)<р(г) , где -6М т/\х(4-р ч 4к *Г)] 3 

т?= ё&у& с]- функция ^ шеет при *“ ><э 

резонансный экстремум в точке 2"= /^/ шириной V э 
а функция на промежутке І/*І+ЫІ} 

практически постоянна , учитывая, что / . Положим 

поэтому ^(г)- </&) и проинтегрируем ч>(Х) : 






Приложение 2 


^ , определенное соотношениями (32)-(34) можно, в 
случав эллипсоидальных образцов с высокосимметричной решеткой, пред- 
СТсШИТЬ в видѳ 


к, = - $ Г - г ^ 

<*■ У V х <а.- 4)І*і-4І 1 

^ 4-<-/ 

* 3 - ^ 7704Г)’ 

а.= ± іні+р х )+ (*-р к )е~*(4 +« 0] - 


(ШО) 

(ПІІ) 

(ПІ2) 

(ПІЗ) 
-], (ПІ4) 




22 


Здесь • а 3 определено при выводе 

(П4). От формы образца зависит только Кр Во всех остальных 
решеточных суммах будни считать, что образец является сферичес¬ 
ким.. Для перехода к ПСС представим в форме 




И 9(<*,->о) ■ 
г ; <к. 


- Л &(4;>0) 


-+ 


(Щ5) 


Для кубических решеток последний член равен нулю. В "низкотемпе¬ 
ратурном" пределе , когда от взаимодействия с ближайши- 

[ зависит только первый член из (ПІ5) , но он не дает 
ж. В остальных членах, а также в К 3 и К 4 
і по угловым и радиальным переменным разделяются 
после замены *1-+ Ъ //3 (і-Зг^Ѳ)І ^ 3 , а легко вычисляется в 
явной форде, что и приводит к результатам, изложенным в основ¬ 


ном тексте. 
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